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1 Abstract 

In this paper we study the higher secant varieties of Grassmann varieties in relation 
to Waring's problem for alternating tensors and to Alexander-Hirschowitz theorem. 
We show how to identify defective higher secant varieties of Grassmannians using a 
probabilistic method involving Terracini's Lemma, and we describe an algorithm which 
can compute, by numerical methods, dim{G{k, n)") for n < 14. Our main result is that, 
except for Grassmannians of lines, if n < 14 and k < (if n = 14 we have studied the 
case fc < 5) there are only the four known defective cases: G(2, 6)^, G{3,7)^, G{3,7)^ 
and G(2,8)''. 

2 Introduction 

Waring's problem for alternating tensors can be expressed in the following form 
(see P) 



Given a vector space V of dimension n + 1 and an alternating tensor 
CO G A'^^^ ^1 what is the least integer s such that lo can be written 
as the sum of s decomposable tensors of the form vi A . . . A ffc+i? 
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This problem is still open and in this paper we will give some evidence for what 
we expect the correct answer to be. 

In order to formulate our result, we will consider a vector space V of dimension 
n + 1 defined over a field IK of characteristic zero, and the Grassmann variety 
G(K''^^, V) = G{k, n), which parametrises the decomposable tensors in the pro- 
jective space of /\^^^ V. As will be explained in the next section, the problem 
translates into finding the dimension of the s-secant variety G{k, nY (see defi- 
nition ESI) • The expected dimension of G(A;, nY is mm{(^J^|) — 1, s(n — k){k + 
1) + s — 1}, otherwise G{k,nY is called defective (see definition VA.l^ . 
It is well known that the Grassmannians of lines G{¥^ ,F^Y defective until 
they fill the ambient space and a list of four defective G{k,nY is given in 0. 
We would like to know if there exist other defective varieties which are still 
unknown. 

Computing the dimension of G{k, nY is quite difficult, even with the aid of a 
symbolic computation package; indeed just after the defective examples of 0, 
the computer's memory reaches its limit with the usual elimination technique 
using Grobner basis. 

The main idea behind this paper is that one can compute dim{G{k,nY) by 
means of a probabilistic method, which consists in studying the span of the 
tangent spaces at s chosen random points. The dimension of this span can be 
computed by numerical methods as the rank of a large matrix, and when this 
dimension coincides with that expected, we can be sure that G{k,nY is not 
defective, indeed with another choice of points the dimension cannot be larger 
because of inequality 

This technique allows us to take the computations further and our main result 
is the following. 

Theorem 2.1. Ifn < 14 and k < ^(ifn = 14 we consider k < 5), G{k,nY 

is defective only for 

^If we choose a basis for V there is a natural 1-1 correspondence between the associated 

bases of /\''^^ V and V, and between the varieties G{k, n) and G(n — fc — 1, n). Thus we 
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• A; = 1 



• {k,n,s) = (2,6,3), S = 1 



• {k,n,s 



) = (3,7,3), 5 = 1 



• {k,n,s) = (3,7,4), S = A 



• {k,n,s 



) = (2,8,4), 5 = 2. 



This theorem is equivalent to the following answer to Waring's problem: 

Corollary 2.2. Given a finite dimension vector space V and an alternating 
tensor lo G A'^^^ if n < 14 and k < (if n = 14 we consider k < 5), then 
u can he written as the sum of s decomposable tensors of the form ui A. . . Aufe+i, 



where s = [^^^^^^i-^ (^+i)] , except for {k, n, s) = (2, 6, 3) , (3, 7, 3) , (3, 7, 4) , (2, 8, 4) 



and k = 1. 

3 Waring's problem and some notations 

Waring 's polynomial problem has attracted considerable attention from geome- 
ters and algebraists throughout its long and absorbing history since it was first 
put forward in 1770. This problem is connected with crucial issues in both re- 
presentation theory and coding theory. 

It poses the following question: if / is a homogeneous polynomial of degree k in 
n variables, what is the least integer s such that / can be written as the sum of 

/cth-powers of s linear forms? 

This formulation of Waring's problem was solved in 1995 by J. Alexander and 
A. Hirschowitz, who produced a formula for finding the integer s; nevertheless 
this formula has four well known exceptions. 

will only consider the variety G{k,n) whore k < 
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Theorem 3.1. 0/ Let char(K) = 0. ^ homogeneous polynomial f € lK[Xo, . . . , X„] 
of degree k can be represented as the sum of s powers of linear forms 

f = L\ + ... + Ll, 

where s = except in the cases where {k,n,s) = (4,2,5) , (4,3,9) , 

(4,4,14) , (3,4,7) and k = 2. 

This challenging result can also be expressed in geometrical terms, as we will 
explain. 

First let us recall some fundamental definitions. 

Let X C be an n-dimensional irreducible projective variety, 

Definition 3.2. The s-secant variety is the closure of the union of all linear 
spaces spanned by s points of X, which is expressed as follows: 

X* = y < xi,...,rc, >. 

Xl,...,Xs£X 

The choice of s points in X gives rise to sn free parameters. In addition, 
s points span a space of projective dimension s — 1 and X'^ will be embed- 
ded in P^. Consequently, we should expect the dimension of X^ to be given by 
min{N, sn + s — 1}: this is called the expected dimension for secant varieties. 

The following estimate on the dimension of X'^ is valid in general: 

dimiX"") < min{iV, sn + s - 1}. (1) 

Prom our viewpoint, the cases where the strict inequality applies are the most 
interesting. 

Definition 3.3. The secant variety X^ is called defective if 

dim{X^) < min{N, sn + s — 1} 
and the quantity 6 = min{N, sn + s — 1} — dim{X^) is its defectiveness. 
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One well-known result that is useful in finding the dimensions of the multi- 
secant varieties is Terracini's lemma. 

Lemma 3.4 ([3]). Let xi,...,Xs ^ X he generic points; let us refer to the 
projectivised tangent spaces to X at these points as T^^X, . . . , T^^X, then 

dim{X') = dim{< T^^X, T^^X >). 

Let us now take a homogeneous polynomial /{Xq, . . . , Xn) of degree k. 
Asking whether / can be written as the sum of powers of degree A; of s linear 
forms Li, . . . ,Ls is the same as asking whether / belongs to the s-secant variety 
of the A:th Veronese embedding of P", which we call Vk^n+i- It is therefore 
important to know the dimension of Vi^^_^_-^, and consequently the cases where 
Vk^n+i is defective for s-secant varieties. 

The result obtained by Alexander and Hirschowitz is extremely useful in our 
case and translates geometrically as follows: 

Theorem 3.5. (l^) The Veronese variety Vk^n+i is defective for s-secant varie- 
ties only in the following cases: 

{k, n, s) = (4, 2, 5) , (4, 3, 9) , (4, 4, 14) , (3, 4, 7) 

and k = 2. 

We have therefore obtained a full classification of defective Veronese varieties. 

In this paper we will analyse the problem of defectiveness with respect to another 
important family of classical varieties, the Grassmannians,^ which are related 
to exterior algebras. 

If A; < dim(y) is a positive integer, we define the Grassmannian G(A;, V) to be 

the variety of projective subspaces of ¥{y) of dimension k. When V = K"^^, 

G{k,V) win be denoted by G{k,n). 

^For the problem of defectiveness of Segre Varieties and its connection with the rank of 

tensors, see [S]. 
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Since dim{G{k,n)) = (k + l){n — k), the expected dimension for the secant 
varieties G{k, nY is: 

min| ^ -|^^ ~ li ^{''T' ~ k){k + 1) + s — l| 
We also have the two following important theoretical results to draw on. 
Theorem 3.6. G{\,ny is defective for s < [|J. 

Theorem 3.7. (;2l) Let k > 2. If s{k + 1) < n + 1 , then G{k,nY has the 
expected dimension. 



4 A probabilistic algorithm and proof of Theoren j2.1l 

To tackle our problem, we initially used the Macaulay 2 computation system 
(see [Z]), which was designed to study problems of algebraic geometry and com- 
mutative algebra, to write an algorithm generating parametric equations for the 
Grassmannians we are studying (see [H]). 

To calculate the dimensions of the multisecant varieties, we favoured a proba- 
bilistic approach involving Terracini's lemma. We took s random points in 
G{k, n) and studied their tangent spaces and the space spanned by these tan- 
gent spaces. If we found the expected dimension, the result was clearly correct, 
but if this revealed defectiveness, more checks needed to be performed. 
Using this approach we constructed an algorithm that turned our problem 
into the calculation of the rank of fairly large matrices with constant coeffi- 
cients; to study dim{G{k, n)^) we needed to know the rank of a matrix of order 
s(l + (/c + l)(n — fc)) X A^. This algorithm enabled us to compute the dimension of 
G{k,ny when n < 11, k <4 and s < 4, at which stage the computer's memory 
was used up. It was therefore clear that symbolic computation was not the best 
tool for this type of task. 

To further proceed with our study, we decided to employ the Matlab software 
system, which is a computation system designed for dealing with numerical com- 
putations involving very large matrices. 
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The new algorithm obtained confirms the vahdity of the probabihstic approach. 
It is based on theoretical observation that the tangent spaces can be computed 
without having to define equations for the Grassmannian. 

Let us take a point P = voA.../\Vk€G = G{k, n); the Plucker coordinates of 
P are all the {k + 1) x {k + 1) minors of the matrix A of order {k + 1) x (n + l), 
which has the vectors vq, . . . ,Vk & V ioi rows. It is easy to check, by the Leibniz 
rule, that the following is true: 

Lemma 4.1. Tp{G) is the projective space associated with 

FAvi A. . .Avk + VQf\V AV2A. . .Avk + • • • + vqA. . .Avfe-iAF = Tq + . . . + Tfe. 

If Ai^j stands for the matrix obtained from A by replacing the zth row by 
the jth row of the identity matrix / of order (n + l) x (n + l), then every Tj is 
parametrised by a matrix Mi of order (n + l) x N whose jth row rrij contains 
the minors of maximum order of the matrix Ai^j. 



( 



AiA 



vo 

Vi-1 

1 ... 

Vi+l 



\ 



( 



vq 

Vi-\ 

... 1 

Vi+\ 



\ 



\ Vk 



Ai,n+1 (2) 



M, 



Our algorithm is described below. 



mi 



Input: positive integers n^ and um- 
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• Repeat on parameters rim < n < um, 1 < A; < [^^J and 2 < s < S = 
r (.+i)(L.)+i a:Dl)to study G^ 

• Define the matrix TA that contains the actual dimensions and the defec- 
tiveness of C. 

• Define the function ed{k, n, s) that calculates the expected dimension of 

and the matrix E of the expected dimensions. 

• Choose s random points Pi, . . . , in G. 

— Take a matrix B of order s{k + 1) x (n + 1) with random rational 
coefficients in the interval [— 

— Extract s submatrices A of order (/c + 1) x (n + 1) from B. 

• Repeat for 1 < /i < s and study Tp^G = Ti + . . . + Tfe+i. 

— Repeat for 1 < i < A; + 1. 

— For every 1 < j < n + 1 calculate the minors (A; + 1) x {k + 1) of 
computed from yl as in ((JJ , and call the row of these minors rrij . 

— Construct the matrix Mi with rows rrij. 

• Parametrise Tp^G + . . . + Tp^G. 

— Concatenate Mi, . . . , M^+i vertically to obtain matrix M. 

• Determine the value of the projective dimension of G*. 

— Calculate the rank of M, then subtract 1. 

— Define row dim of TA of actual dimensions and row dif = E — dim 
of defectiveness. 

• Output: matrix TA. 

This is the text of the algorithm. 
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L=100 

nin=3 

nM=14 

f=irLline( 'floor ( (n-l)/2) ') 
F=[] 

for u=iim:nM 
F=[F 2*f (u)] 
end 

N=inline ( 'factorial (n+1) / (factorial (k+1) *f actorial (n-k) ) -1 ' ) 
S=inline ( ' ceil ( (N+1) / ( (n-k) * (k+1) +1) ) ' ) 
Smax=S(N(f (nM) ,nM) ,f (nM) ,nM) 
TA=[] 

for n=nm:nM 
T=zeros(2*f (n) ,Smax+l) 
for k=l:f(n) 
E=[] 

T(2*k-l,l:2) = [k n] 
I=eye(n+1) 

v=nchoosek(l :n+l ,k+l) 

l=size(v, 1) 

dim=[] 

for s=2:S(N(k,n) ,k,n) 

M=[] 

ed=inline ( ' min (N , s* (n-k) * (k+1 ) +s- 1 ) ' ) 
E=[E ed(N(k,n) ,k,n,s)] 
B=rand(k+1, (n+l)*s) 
B=(B-0.5)*2*L 
for h=l:s 

A=B(: , (h-l)*(n+l)+l:h*(n+l)) 
for i=l:k+l 
r=A(i,:) 

for j=i:n+i 

A(i,:)=I(j,:) 

m=[] 

for w=l:l 

D(w)=det(A(: ,v(w, :))) 
m=[m D(w)] 
end 
M=[M;m] 
end 
A(i,:)=r 
end 
end 

dim=[dim rank(M)-l] 
dif=E-dim 
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T(2*k-1 , s+l)=dim(l , s-1) 
T(2*k,s+l)=<iif (l,s-l) 
end 
end 

TA=[TA;T] 
end 

TA 

Us = S = r (fc+i)(n-fc)+i (fc+i)l G{k, ny is not defective, then the variety 
fills the ambient space = p(fc+i)^^; if it is defective, then we find dif > 
and it can happen that dim{G{k,n)^) < N, so that we will have to calculate 
dim{G{k,ny) for s> S. 

Using this algorithm, at the stage (n, A;) = (6, 14) the computer's memory was 
used up. Our results are summarised in the following tables. 
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These results confirm what is known about the defectiveness of Grassman- 
nians. Except for Grassmannians of hues, we have identified four defective 
varieties. We have therefore proved theorem 12.11 
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